1 Introduction

  1. TMLE is a general algorithm for the construction of double robust, semiparametric, efficient substitution estimators. TMLE allows for data-adaptive estimation while obtaining valid statistical inference.
  2. Although TMLE implemtation uses the G-computation estimand (G-formula). Briefly, the TMLE algorithm uses information in the estimated exposure mechanism P(A|W) to update the initial estimator of the conditional mean E_0(Y|A,W).
  3. The targeted estimates are then substituted into the parameter mapping. The updating step achieves a targeted bias reduction for the parameter of interest \(\psi\)(P0) (the true target parameter) and serves to solve the efficient score equation. As a result, TMLE is a double robust estimator.
  4. TMLE it will be consistent for \(\psi\)(P0) is either the conditional expectation E_0(Y|A,W) or the exposure mechanism P_0(A|W) are estimated consistently. When both functions are consistently estimated, the TMLE will be efficient in that it achieves the lowest asymptotic variance among a large class of estimators. These asymptotic properties typically translate into lower bias and variance in finite samples.(Bühlmann et al., 2016)
  5. The advantages of TMLE have been repeatedly demonstrated in both simulation studies and applied analyses.(Laan and Rose, 2011)
  6. The procedure is available with standard software such as the tmle package in R (Gruber and Laan, 2011).

2 The G-Formula

  1. \(\Psi(P0)\,=\,\sum_{w}\,\left[\sum_{y}\,P(Y=y\mid A=1,W=w)-\,\sum_{y}\,P(Y = y\mid A=0,W=w)\right]P(W=w)\)
    where
    \(P(Y = y \mid A = a, W = w)\,=\,\frac{P(W = w, A = a, Y = y)}{\sum_{y}\,P(W = w, A = a, Y = y)}\)
    is the conditional probability distribution of Y = y, given A = a, W = w and,
    \(P(W = w)\,=\,\sum_{y,a}\,P(W = w, A = a, Y = y)\)
  2. Using classical regression methods to control confounding requires making the assumption that the effect measure is constant across levels of confounders included in the model.
  3. Alternatively, standardization allows us to obtain an unconfounded summary effect measure without requiring this assumption.The G-formula is a generalization of standardization(Greenland and Robins, 1986)

3 TMLE flow chart

Source: (???)

4 Implementation

4.1 First step

5 Thank you

Thank you for participating in this tutorial.

If you have updates you would like to make to the lesson, please send me a pull request.
Alternatively, if you have any questions, please e-mail me.
Miguel Angel Luque Fernandez
E-mail: miguel-angel.luque at lshtm.ac.uk
Twitter @WATZILEI

6 Session Info

devtools::session_info()

7 References

Bühlmann P, Drineas P, Laan M van der, Kane M. (2016). Handbook of big data. CRC Press.

Greenland S, Robins JM. (1986). Identifiability, exchangeability, and epidemiological confounding. International journal of epidemiology 15: 413–419.

Gruber S, Laan M van der. (2011). Tmle: An r package for targeted maximum likelihood estimation. UC Berkeley Division of Biostatistics Working Paper Series.

Laan M van der, Rose S. (2011). Targeted learning: Causal inference for observational and experimental data. Springer Series in Statistics.

LS0tCnRpdGxlOiAiVE1MRSBzdGVwIGJ5IHN0ZXAiCmF1dGhvcjogJ0J5OiBNaWd1ZWwgQW5nZWwgTHVxdWUgRmVybmFuZGV6LCBtaWd1ZWwtYW5nZWwubHVxdWVAbHNodG0uYWMudWsnCmRhdGU6ICJPY3RvYmVyIDE1dGgsIDIwMTYiCm91dHB1dDogIAogIGh0bWxfbm90ZWJvb2s6CiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIGhpZ2hsaWdodDogZGVmYXVsdAogICAgI2tlZXBfbWQ6IHllcwogICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMKICAgIHRoZW1lOiBqb3VybmFsCiAgICB0b2M6IHllcwogICAgdG9jX2Zsb2F0OgogICAgICBjb2xsYXBzZWQ6IG5vCiAgICAgIHNtb290aF9zY3JvbGw6IHllcwogICAgICB0b2NfZGVwdGg6IDMKY3NsOiByZWZlcmVuY2VzL2lzbWUuY3NsCmJpYmxpb2dyYXBoeTogcmVmZXJlbmNlcy9iaWJsaW9ncmFwaHkuYmliCmZvbnQtaW1wb3J0OiBodHRwOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1SaXNxdWUKZm9udC1mYW1pbHk6ICdSaXNxdWUnCi0tLQoKPCEtLUJFR0lOOiAgU2V0IHRoZSBnbG9iYWwgb3B0aW9ucyBhbmQgbG9hZCBwYWNrYWdlcy0tPgpgYGB7ciBzZXQtZ2xvYmFsLW9wdGlvbnMsIGVjaG8gPSBGQUxTRX0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGV2YWwgPSBUUlVFLCAKICAgICAgICAgICAgICAgICAgICAgIGVjaG8gPSBUUlVFLCAKICAgICAgICAgICAgICAgICAgICAgIGNhY2hlID0gRkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICBpbmNsdWRlID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgIGNvbGxhcHNlID0gRkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICBkZXBlbmRzb24gPSBOVUxMLAogICAgICAgICAgICAgICAgICAgICAgZW5naW5lID0gIlIiLCAjIENodW5rcyB3aWxsIGFsd2F5cyBoYXZlIFIgY29kZSwgdW5sZXNzIG5vdGVkCiAgICAgICAgICAgICAgICAgICAgICBlcnJvciA9IFRSVUUsCiAgICAgICAgICAgICAgICAgICAgICBmaWcucGF0aD0iRmlndXJlcy8iLCAgIyBTZXQgdGhlIGZpZ3VyZSBvcHRpb25zCiAgICAgICAgICAgICAgICAgICAgICBmaWcuYWxpZ24gPSAiY2VudGVyIiwgCiAgICAgICAgICAgICAgICAgICAgICBmaWcud2lkdGggPSA3LAogICAgICAgICAgICAgICAgICAgICAgZmlnLmhlaWdodCA9IDcpCmBgYAoKI0ludHJvZHVjdGlvbgoxLiAqKlRNTEUqKiBpcyBhIGdlbmVyYWwgYWxnb3JpdGhtIGZvciB0aGUgY29uc3RydWN0aW9uIG9mIGRvdWJsZSByb2J1c3QsIHNlbWlwYXJhbWV0cmljLCBlZmZpY2llbnQgc3Vic3RpdHV0aW9uIGVzdGltYXRvcnMuIFRNTEUgYWxsb3dzIGZvciBkYXRhLWFkYXB0aXZlIGVzdGltYXRpb24gd2hpbGUgb2J0YWluaW5nIHZhbGlkIHN0YXRpc3RpY2FsIGluZmVyZW5jZS4gCjIuIEFsdGhvdWdoICoqVE1MRSoqIGltcGxlbXRhdGlvbiB1c2VzIHRoZSBHLWNvbXB1dGF0aW9uIGVzdGltYW5kIChHLWZvcm11bGEpLiBCcmllZmx5LCB0aGUgVE1MRSBhbGdvcml0aG0gdXNlcyBpbmZvcm1hdGlvbiBpbiB0aGUgZXN0aW1hdGVkIGV4cG9zdXJlIG1lY2hhbmlzbSBQKEF8VykgdG8gdXBkYXRlIHRoZSBpbml0aWFsIGVzdGltYXRvciBvZiB0aGUgY29uZGl0aW9uYWwgbWVhbiBFXzAoWXxBLFcpLiAKMy4gVGhlIHRhcmdldGVkIGVzdGltYXRlcyBhcmUgdGhlbiBzdWJzdGl0dXRlZCBpbnRvIHRoZSBwYXJhbWV0ZXIgbWFwcGluZy4gVGhlIHVwZGF0aW5nIHN0ZXAgYWNoaWV2ZXMgYSB0YXJnZXRlZCBiaWFzIHJlZHVjdGlvbiBmb3IgdGhlIHBhcmFtZXRlciBvZiBpbnRlcmVzdCAkXHBzaSQoUDApICh0aGUgdHJ1ZSB0YXJnZXQgcGFyYW1ldGVyKSBhbmQgc2VydmVzIHRvIHNvbHZlIHRoZSBlZmZpY2llbnQgc2NvcmUgZXF1YXRpb24uIEFzIGEgcmVzdWx0LCBUTUxFIGlzIGEgKipkb3VibGUgcm9idXN0IGVzdGltYXRvcioqLgo0LiAqKlRNTEUqKiBpdCB3aWxsIGJlIGNvbnNpc3RlbnQgZm9yICRccHNpJChQMCkgaXMgZWl0aGVyIHRoZSBjb25kaXRpb25hbCBleHBlY3RhdGlvbiBFXzAoWXxBLFcpIG9yIHRoZSBleHBvc3VyZSBtZWNoYW5pc20gUF8wKEF8VykgYXJlIGVzdGltYXRlZCBjb25zaXN0ZW50bHkuIFdoZW4gYm90aCBmdW5jdGlvbnMgYXJlIGNvbnNpc3RlbnRseSBlc3RpbWF0ZWQsIHRoZSAqKlRNTEUqKiB3aWxsIGJlIGVmZmljaWVudCBpbiB0aGF0IGl0IGFjaGlldmVzIHRoZSBsb3dlc3QgYXN5bXB0b3RpYyB2YXJpYW5jZSBhbW9uZyBhIGxhcmdlIGNsYXNzIG9mIGVzdGltYXRvcnMuIFRoZXNlIGFzeW1wdG90aWMgcHJvcGVydGllcyB0eXBpY2FsbHkgdHJhbnNsYXRlIGludG8gbG93ZXIgYmlhcyBhbmQgdmFyaWFuY2UgaW4gZmluaXRlIHNhbXBsZXMuW0BidWgyMDE2XQo1LiBUaGUgYWR2YW50YWdlcyBvZiBUTUxFIGhhdmUgYmVlbiByZXBlYXRlZGx5IGRlbW9uc3RyYXRlZCBpbiBib3RoIHNpbXVsYXRpb24gc3R1ZGllcyBhbmQgYXBwbGllZCBhbmFseXNlcy5bQHZhbjIwMTFdCjYuIFRoZSBwcm9jZWR1cmUgaXMgYXZhaWxhYmxlIHdpdGggc3RhbmRhcmQgc29mdHdhcmUgc3VjaCBhcyB0aGUgKip0bWxlKiogcGFja2FnZSBpbiBSIFtAZ3J1YmVyMjAxMV0uCgojVGhlIEctRm9ybXVsYQoxLiAkXFBzaShQMClcLD1cLFxzdW1fe3d9XCxcbGVmdFtcc3VtX3t5fVwsUChZPXlcbWlkIEE9MSxXPXcpLVwsXHN1bV97eX1cLFAoWSA9IHlcbWlkIEE9MCxXPXcpXHJpZ2h0XVAoVz13KSQgIAp3aGVyZSAgCiRQKFkgPSB5IFxtaWQgQSA9IGEsIFcgPSB3KVwsPVwsXGZyYWN7UChXID0gdywgQSA9IGEsIFkgPSB5KX17XHN1bV97eX1cLFAoVyA9IHcsIEEgPSBhLCBZID0geSl9JCAgCmlzIHRoZSBjb25kaXRpb25hbCBwcm9iYWJpbGl0eSBkaXN0cmlidXRpb24gb2YgWSA9IHksIGdpdmVuIEEgPSBhLCBXID0gdyBhbmQsICAKJFAoVyA9IHcpXCw9XCxcc3VtX3t5LGF9XCxQKFcgPSB3LCBBID0gYSwgWSA9IHkpJCAgCjIuIFVzaW5nIGNsYXNzaWNhbCByZWdyZXNzaW9uIG1ldGhvZHMgdG8gY29udHJvbCBjb25mb3VuZGluZyByZXF1aXJlcyBtYWtpbmcgdGhlIGFzc3VtcHRpb24gdGhhdCB0aGUgZWZmZWN0IG1lYXN1cmUgaXMgY29uc3RhbnQgYWNyb3NzIGxldmVscyBvZiBjb25mb3VuZGVycyBpbmNsdWRlZCBpbiB0aGUgbW9kZWwuCjMuIEFsdGVybmF0aXZlbHksICoqc3RhbmRhcmRpemF0aW9uKiogYWxsb3dzIHVzIHRvIG9idGFpbiBhbiB1bmNvbmZvdW5kZWQgc3VtbWFyeSBlZmZlY3QgbWVhc3VyZSB3aXRob3V0IHJlcXVpcmluZyB0aGlzIGFzc3VtcHRpb24uVGhlICoqRy1mb3JtdWxhKiogaXMgYSAqZ2VuZXJhbGl6YXRpb24gb2Ygc3RhbmRhcmRpemF0aW9uKltAcm9iaW5zMTk4Nl0KCiNUTUxFIGZsb3cgY2hhcnQgCiFbXShGaWd1cmVzL3RtbGUucG5nKQpTb3VyY2U6IFtAdmFuMjAxNl0KCiNJbXBsZW1lbnRhdGlvbiAKIyNGaXJzdCBzdGVwCgoKIyBUaGFuayB5b3UgIApUaGFuayB5b3UgZm9yIHBhcnRpY2lwYXRpbmcgaW4gdGhpcyB0dXRvcmlhbC4gIAoKSWYgeW91IGhhdmUgdXBkYXRlcyB5b3Ugd291bGQgbGlrZSB0byBtYWtlIHRvIHRoZSBsZXNzb24sIHBsZWFzZSBzZW5kIDxhIGhyZWY9Imh0dHBzOi8vZ2l0aHViLmNvbS9taWdhcmlhbmUvTUFMRiIgdGFyZ2V0PSJfYmxhbmsiPm1lPC9hPiBhIHB1bGwgcmVxdWVzdC4gIApBbHRlcm5hdGl2ZWx5LCBpZiB5b3UgaGF2ZSBhbnkgcXVlc3Rpb25zLCBwbGVhc2UgZS1tYWlsIG1lLiAgCioqTWlndWVsIEFuZ2VsIEx1cXVlIEZlcm5hbmRleioqICAKKipFLW1haWw6KiogKm1pZ3VlbC1hbmdlbC5sdXF1ZSBhdCBsc2h0bS5hYy51ayogIAoqKlR3aXR0ZXIqKiBgQFdBVFpJTEVJYCAgCgojIFNlc3Npb24gSW5mbyAKYGBge3Igc2Vzc2lvbi1pbmZvLCByZXN1bHRzID0nbWFya3VwJ30KZGV2dG9vbHM6OnNlc3Npb25faW5mbygpCmBgYAojIFJlZmVyZW5jZXMgCg==